On the integrability of the square-triangle random tiling model 



On 



O 



Jan de Gier, Bernard Nienhuis 
Instituut voor Theoretische Fysica, Universiteit van Amsterdam, Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 

(February 9, 2008) 

It is shown that the square-triangle random tiling model is equivalent to an asymmetric limit of 
the three coloring model on the honeycomb lattice. The latter model is known to be the 0(n) model 
at T = and corresponds to the integrable model connected to the affine A% Lie algebra. Thus it is 
shown that the weights of the square-triangle random tiling satisfy the Yang-Baxter equation, albeit 
in a singular limit of a more general model. The three coloring model for general vertex weights is 
solved by algebraic Bethe Ansatz. 
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I. INTRODUCTION 



> 

Random tiling models have gained renewed interest in the past years by the discovery of quasi-crystals p|. They 
provide an example of the entropic stability of structures whose diffraction pattern has a rotational symmetry which 
is incompatible with periodicity. As such they offer an explanation of the existence of quasi-crystalline alloys ||. 
In analogy with the diamond covering used to describe the ground state configurations of the triangular Ising anti- 
I ■ ferromagnet ||], random tilings can be described by a domain wall structure. The difficulty is that there is more 
than one type of domain wall, in contrast to the diamond covering. This fact makes it much more difficult to solve 
the model by coordinate Bethe Ansatz, as was tried for the octagonal square-rhombus tiling Q. Widom Q however 
succeeded in diagonalizing the transfer matrix using Bethe Ansatz for the 12-fold square-triangle random tiling model. 
Shortly after that Kalugin (|] was able to find a closed expression for the entropy as a function of the domain wall 
densities in part of the phase diagram. More recently a solution of an octagonal random tiling has also been found 

Being solvable by coordinate Bethe Ansatz, it is natural to look for solutions of the Yang-Baxter equation RBI for 
these models. This would provide an answer to why these models are integrable and a canonical way to diagonalizc 
the transfer matrix via the algebraic Bethe Ansatz ||. The matrix of Boltzmann weights j§] however, is not invertible, 
which is necessary to obtain commutativity of the transfer matrix. In this paper we show that the square-triangle 
random tiling model is a singular limit case of a more general model obeying the Yang-Baxter equation. 



II. THE MODEL 



We shall consider a vertex model on the square lattice, whose Boltzmann weights are denoted by 

. p. : 

W((J,,a;/3,v\u) = M — 



(1) 



Each edge of the lattice can be in one of three different states, 1,2 and 3. The partition function of the model is given 
by 

Z = E U w ®> ( 2 ) 

config. i 

where we sum over all configurations which are weighted by the product over the vertices i of their local Boltzmann 
weights W(i). The explicit form of the weights can be found in Table |[ The model can be written in terms of the 
weights Wq([j,, a; [3, v) associated with the affine Lie algebra which can be found in |l|. 

W(n,a;P,v\u) = O^,W {n' ,a; 0,v'\u)O v , v , (3) 
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where O = diagjxi, x%, X3}. The weights ^(a, v, P, v) satisfy the star-triangle or Yang-Baxter equation (YBE) 8, 

Y, W (H, v; z/» - u)W (n", a; 7, M»Wb(f", 7; P, f» = 

7,1*", v " 

Y W (p,a- n ,v''\u)W ( l x, T ,p, l j,''\v)W (^,v''; V ', l j,'\v-u). (4) 

7,1*","" 

Equation (|j) can be written more elegantly by defining the operators L(u) 

( L ( u )^) a i3 = W (fi, a; (3, v\u), 

L a (u) : V a <g> C 3 -^V a ® C 3 , V a ~ C 3 . 

The auxiliary label a is introduced for later convenience. We shall omit it if there is no confusion on which space L{u) 
acts. In this language the vertex states a (a = 1, 2, 3) are represented by the standard basis e a of C 3 . V a — Vb — C 3 are 
so-called auxiliary spaces, corresponding to the horizontal edges of a vertex. The transfer matrix T = 52 M=1 



on a lattice of horizontal size N can be written in terms of the local operators L a (u) (IJ) as 

/ n ^ 

Ta(u)= (n°a^>) 

N (6) 

T„(«) : K <g> W -» V a <g> W, W - (g) C 3 , 
T(u) = Tr a T a (w). 

The trace is taken only over the auxiliary space matrix structure. The operators L 3 a (u) act as L a (u) on the jth factor 
in 7i and as the identity on all other factors, 

j — 1 times N— j times 

L J (u)^ =1 <g> • • • (8 1® L(u) M „!g) 7 <g> • • • <8> I, (7) 
where I is the identity on C 3 . The partition function (||) on a lattice of size N x M can then be written as 

Z = Tr w 7» M . (8) 

Furthermore we define the i?-matrix as 

(R(v - u)v, u ) a0 = W ((i, «; P, v\v -u), 

R ab (v -u): K®14-»K® H- <J1 

Here too, the Roman labels a and 6 indicate on which auxiliary space R is acting. Greek labels will be used to indicate 
matrix elements. The YBE (|]) can be written as 

(R(v - u)^») vv n [L(u) M // M 'L(u) l ,« !/ /] = 
[L(u) (10) 

where summation over repeated indices is understood. Each element L{u)a V of £(u) is an operator acting on C 3 . We 
shall regard R(u — v) as a 9 x 9 matrix which is given explicitly in Table |. In a compact notation the YBE ( |l0| ) can 
be written as an operator equation on the tensor product V a ® Vb ® C 3 . It then assumes the guise 

Rab(v - u) ■ [L a (v) <g> L 6 (tt)] = [L b (u) ® • i? Qb (u - it). (11) 

From ( [TT| ) we obtain the YBE for the matrix T(u) as defined in (JsJ) . 

i? ab (« - u) ■ [T a (v) <g. T 6 («)] = [T 5 (u) (8 T a («)] • R ab {v - u). (12) 

From the expression (|J) for the partition sum one sees that the leading term is given by the largest eigenvalue of 
T. For the sake of completeness we give the explicit diagonalization of the transfer matrix using the algebraic Bethe 
Ansatz in appendix |A| This is just a special case of the trigonometric case of Kulish and Reshetikhin p3[ . 
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W(fj,, (i; /tx, fx\u) — x 2 ^ sinh(u + A), 

W(n, u; n, v\u) = x^x„c- u "S"^-") sin h(A), 

W{1, 2; 2, = x\yl 2 sinh(u), W{1, 3; 3, = x\y^ 2 sinh(u), 

W{2, 3; 3, 2\u) = x\y^ 2 sinh(u), W(2, 1; 1, 2|u) = a&/§ sinh(n), 

W{3, 1; 1, 3|u) = x\y\ sinh(u), W(3, 2; 2, 3|u) = i§j/? sinh(u), 



/S2 
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TABLE I. Boltzmann weights W and i?-matrix corresponding to the weights Wo. Here we use the abbreviations so = sinh(A), 
si = sinh(u), S2 = sinh(u + A). 



III. HONEYCOMB LATTICE 



As was already shown by Reshetikhin [^|, at a special value of the spectral parameter u, the model defined by 
the weights Wq factorizes on the honeycomb lattice. In this section we rederive this result for the model in Table |. 
Consider the operators 



OL{u)Ot : C 3 (g>V ^ V< 



(13) 



where r is the permutation operator in C 3 ® C 3 (recall that V ~ C 3 ), r(e Q ® ep) — ep ® e a . The eigenvectors of 
OL(u)Ot are given by 



Vl + e" 



2\ 



Vl + C 



(xpy~ t ep®e a -y 1 1 x a e A e Q ®e (3 ), 
=(y^Xpep ® e Q + y~ 1 x a e x e a <g> e^), 



/(77) — e 7 <8> e 7 , 

with (a, (3, 7) a cyclic permutation of (1,2,3). These vectors satisfy the eigenvalue equations 

OL(u)Ore 7 = —x a xp sinh(it — A)e 7 , 
OL(u)Orf^ v ) = x^x v sinh(w + A)/( M „). 



Thus at M = — A the operator Lt becomes a projector. 










sinh(2A)E n 



/' 



(14) 



(15) 



FIG. 1. Factorization of the weights 
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Introducing the dual vectors e* with the properties 

e* • e 7 , = XaXpSj*/', (16) 

we can write OLOt at u = —A as 

3 

OL(-X)Ot = sinh(2A) ^ e 7 e 7 . (17) 

7=1 

Give the y- vertices of the honeycomb lattice the weights e* • (ep ® e„) and the A- vertices the weights (e^ ® e Q ) • e 7 . 
Graphically, the matrix element (£(— A)r /i/ g) Qly , corresponding to the vertex weight W(/x, a; /3, i/| — A), can be written 
as a sum over products of two vertices of the honeycomb lattice, see Figure 0. 

It follows that the model on the honeycomb lattice with vertex weights given in Figure || (where (123) — > (ABC)) 
has the same partition function as the model in Table | at u = —A on the square lattice. More precisely, 

S mh(X)- NM Z s N Q M (u = -A) = Z?C M (\), (18) 



^A B^ .A ^B A^ 









C A B C A B 

2V3C' X/2 -x 3 y lC - x/2 x 1 y- 1 c x/2 xxy' 1 ^' 2 x 2 y- 1 c x/2 

ab -1 bc^ 1 ca —1 ba~ 1 cb~ 1 



C A B C A B 









B A C B A' A B B C C A 



-1 — A/2 -1 — A/2 -1 -A/2 A/2 

-xii/ 2 o ' —xty 3 e ' -X23/1 e x 3 y 2 e ' 



cb- 1 



FIG. 2. Vertex configurations on the honeycomb lattice. On the first line their weights are given and on the second line the 
gauge transformations. 

The partition function does not change if we apply a gauge transformation to the weights. If we multiply the 
weights by the factors shown in Figure || and choose b = — e _A / 3 a and c = e~ 2A / 3 a and set Xi = j/j = 1, we obtain the 
partition sum of the fully packed loop (FPL) model on the honeycomb lattice ]l4| , 

z s N Q M (u = -\) = shMX) NM Y, nNid) - ( 19 ) 

c 

The sum runs over all dense loop coverings C of the honeycomb lattice. N(C) is the number of loops in the covering 
C and n = 2 cosh(A) is the loop fugacity. 

The edge states may be interpreted as differences modulo three (going clockwise around each vertex) between three- 
state Potts variables on the vertices of the triangular lattice. If the states A, B and C are interpreted as differences 0,1 
and 2, the corresponding Potts model allows only configurations on the triangle in which two variables are equal and 
the third is different. The fully ferromagnetic and the completely anti- ferromagnetic arrangements are then excluded. 
This model has competing ferromagnetic two-spin and anti- ferromagnetic three-spin interactions. 

IV. SQUARE-TRIANGLE TILING 

Take the dual of the honeycomb lattice and associate to each edge of this triangular lattice the corresponding state 
variable of the honeycomb lattice. Relabel the states of the horizontal axis by A, B, C — ► 0, +, — , on the ascending 
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diagonal by A, B, C — > +, — , and on the descending diagonal by A, B, C — > — , 0, +. The partition sum ( |l8| ) is 
then equal to the partition sum of a model on the triangular lattice with face configurations given in Figure |3[ 








+ 






0/\0 







o\ /o 



FIG. 3. Face configurations on the triangular lattice corresponding to the vertex configurations in Figure g 

Now let x^ 1 — X2 = x$ = x 1 / 2 and y\ — y2 — Vz = x^ 1 ^ 2 . If x = the faces with a on all three edges vanish. In 
that case the states + and — can be regarded as rotation angles of edges w.r.t. a fixed triangular lattice. Take these 
angles ±7r/12 and wipe out every edge with state 0. In this way the model maps onto the square triangle random 
tiling model. The state corr espo nds to the diagonal of a square. 

The eigenvalue expression ( A2§| ) in the limit u = — A becomes 



A(-A) = sinh(-A)" ft fl 8i ° h(u %+ 2A) 

k=i smh K + A ) 1=1 smh(uj + A) 

ri2 . (2)\ 

+ x»' sinh(-A)" J] S } ■ (20) 

;=1 sinh(it; ; + A) 



The Bethe Ansatz equations ( A26 ) and ( A27 ) for the two sets of momenta become 



k 

and 



j* sinh(,f -«W) - sinh(,f - + A) = ^ ^ 
x sinh(tt^ — + A) smh(uy - u| 2 ' - A) 



sinh^f + A) ^ _ j;jV+ra2 -^ sinh(4 1) -4 1) + A) 
, sinh(wj 1) ) / k=1 sinh(u' 1) - - A) 

sinh(uj 1 ' — up' — A) 

. w (1) (2)7 

1=1 sinh(u) ; - u\ J ) 



- sinn ti — u, ■ — a 1 



It must be noted that these equations can also be derived from a coordinate Bethe Ansatz. For the FPL model this 
was done by Baxter whose method can be generalized slightly to obtain the equations ([H]) and (|2^). For our 
purposes, taking the limit as — *■ 0, it is more convenient to use Baxter's variables. This can be accomplished by making 
the following substitutions 

_ xsmh(uf ) ) ^ xsmhju^) 

Sj ~ sinh(uf ) + A) ' M ~ sinh(u| I 2) + A) ' 

On taking the limit x — > we arrive at the following expressions for the Bethe Ansatz equations and the eigenvalue 
of the transfer matrix of the square-triangle tiling in its triangular lattice representation 



b 3 



fe=l 

/ "1 



A= [i-r + l[s k )l[t u , (25) 

\ k=l j v=l 

These equations can be solved analogously to the original solution of the square-triangle random tiling by Kalugin 



5 



V. CONCLUSION 



We have made a connection between the recently solved square-triangle random tiling model and a known solvable 
lattice model. It follows from Table | and the substitution = x" 1 / 2 that the i?-matrix for the square-triangle tiling 
(x = 0) is either singular or contains infinite elements. As a result the transfer matrix of the square-triangle tiling 
model is a limit of a family of commuting transfer matrices, but is itself not a member of such family. For any finite 
x though the i?-matrix is invertible, which implies integrability via the Yang-Baxter equation for this more general 
model. The square-triangle random tiling thus is a singular limit of a model which is integrable in the usual sense. To 
obtain the weights of the square-triangle model one has to take the limit u — » — A first and then take x — > 0. These 
two limits do not commute. 

One final point can be made about the robustness of integrability. The square-triangle tiling has been solved in 
three different ways. One is that of Widom and Kalugin (5j|(|, the second can be found in |16) and the third one 
in this paper. All these methods only differ in their choice of representation which of course should not influence 
integrability. 
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APPENDIX A: DIAGONALIZATION OF THE TRANSFER MATRIX 



We write the transfer matrix T{u) as an operator on (V^ © V^) tgiH^V ®H: 



A(u) B(u) 
C{u) D{u) 

The entries of T(u) act on the following spaces: 

A(u) : F (0) ®W-> V(°) <g> H, B{u) : © H 

C{u): ®Ti->V^ ®H, D(u): K« © ft -> © ft. 



(Al) 



(A2) 



~ C and ~ C 2 are the subspaces of V corresponding to the natural embedding of C 2 in C 3 . The R- matrix 
has the following form on the standard basis of (Va ^ © V^ 1 "*) © © V b ^)i 



R a b(u) 



/sinh(u + A) \ 

sinh(u)?7 6 _1 e M sinh(A)I (1 X ) 

e~ u sinh(A)/(°^ 1 ) sinh(w)[/ a 
V R$(u)J 



(A3) 



Here, F a v 



(o) 



V, 



(i) ' 



(0)(1) 



Va (1) © K 



(0) / 



(1)(0) 



Va © ar e canonical isomorphisms and 



r(l) 



2/1 

o v\ 



(A4) 



With the notation U a we denote the operator that acts as U in the space and trivial everywhere else. The 

Oft ■ Va ^ ® Vfr ' II > /, I"" im OfUIJA .">l 1. I LV l I U V <!.- ua I UII / l~UJ 



reduced matrix r[)} : © V^ 1 



/sinh(u + A) 

y~ 2 sinh(u) e" sinh(A) 

e~ u sinh(A) y\ sinh(u) 

V sinh(u + A) / 



(A5) 



From the relation (|12|) all commutation relations between the operators A, B, C and D can be obtained. In the sequel 
we shall only need the following three of them, 



A b (u) © B a (v) = smh(v - it) 1 (sinh(v -it + X)B a (v) © A b (u) 
-e v - u smh^)B b {u)(g,A a {v)I (1 ^ U~\ 

D a (v) © B b {u) = sinh(i; - u)- 1 ^ 1 (B b (u) © D a (v)R$(v - u) 

-e u ~ v sinh(A)/ (0)(1) S a (w) © D b (u) 



B a {u) © B b (v) = smh(v -u+X) L B b (v)<gi B a (u)R$(v - u) 



(A6) 

(A7) 
(A8) 



The first step in diagonalizing T is to construct candidates for its eigenvectors. This construction will be outlined 
below. One first defines a 'pseudo-vacuum ' on which T(u) is upper trigonal, 



N 



(A9) 







It then follows that 



T a (u)F(°> = 



'x\ N sinh(u + A) W F(°) * * 
[xlyl) N smk{u) N * 
" (sii^)^ sinh(u) Ar ^( )/ 



(A10) 



7 



therefore is an eigenvector of T(u). To obtain more eigenvectors we make the following Ansatz 



F = B ni (<>) 



(All) 



The vectors F^'™ 1 will be found later. Schematically, we can represent the action of the transfer matrix on F as an 
(rii |l)xJV lattice, see Figure f|. The state F^ is represented by N vertical edges and the action of each of the 
Bi{u^) by a horiz onta l line. The action of D(u) on the resulting vector F is given by the upper horizontal line. The 
commutation rule ( A7) can now be represented graphically by shi fting the upper line downwards using the YBE (Jlj" 
In a similar fashion we can get a graphical representation of (A6). 



i?(0) 



D : 


y(l) 


B : 


yW 


B : 


yiD 



(i) 



(0) 



(0) 



FIG. 4. Commutation rule (A7) for m = 2. 



The factors that arise after commutation are given in Figure ^ by the external vertices. For example, the vertex on 
the left of the second diagram of Figure |i| co rresponds to (u—u^) and the one on the right to sinh(u — )~ 1 U^' 1 . 
It follows from the relations ( |A6| ) and ( |A7| ) that 



A{u)F = x\ N sinh(u + \) N Yl 



smh(ui — u + A) 



k=i sinh(4 
B ni (u$) <8 ... ® Bxiu^p^ 1 ® • • • ® C/f 1 ^'" 1 
'unwanted terms'. 



(A12) 



= sinh(u - uj^)-^^) ® . ..® SiCu^) x 

k=l 

T a 1) (u;{u ( k 1) })F {1) ' ni + 'unwanted terms'. 
The reduced transfer matrix is given by 

ri 1J (u; {4 1} }) = U-^D a (u)R^(u - u[ iy ) . . . (u - u&>), 



(A13) 



(A14) 



The 'unwanted terms' are similar to (A12) and (A13) but now one of the Bj has u instead of Uj as its argument. 
Provided that the 'unwanted terms' vanish, the vector F will thus be an eigenvector of T(u) = A(u) + Tr a D a (u) with 
eigenvalue A(u), given by 
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A( U) = x r + a) w ft ^ { i\- u \ x) tiu) 

fe=l smh w - u ) 

ni 

+ Y[smh(u-u i k 1) )- 1 A^(u). 

k=l 



(A15) 



J/" 1 ® • • • ® I7f 1 and T« («; {ujj. J }) are diagonal with eigenvalues /x(?7) and AW(tt) respectively. For the 'unwanted 



terms', it is easily seen that the terms in which B nx (v!n}) is replaced by B ni (u) are (up to a common multiplicative 
factor) precisely of the form ( A12| ) and (A13) with u and Un} interchanged and the factor with k — n\ omitted. Using 
the commutation rule for B(u)B{v) it can be shown that the same also holds for the terms with Bj(u^) replaced by 
Bj{u). The 'unwanted terms' will therefore cancel if A.(v,j) = 0, or 

m 

xl N fi(U) Binhfu^ + X) N Y[ sinh(uf 5 - uj^ — A) = -A^uf 5 ) (A16) 

fc=i 

The problem of finding the eigenvectors F^'™ 1 of is completely analogous to the construction above. It follows 
from the commutation rule 



R§j(v - u)D a (v) ® D 6 (u) = D 6 (u) ® A^it^u - «) 



(A17) 



and the fact that R ab has the same form as R a b that obeys the YBE 

Tl 1] (u; {u^}) ® T« (v; {u^})} ■ R§ (v-u). 



Writing TW(u;{uP}) as 



^w^)=(;|;!2|:!) 



(A18) 



(A19) 



since R^(u) has the same structure as R(u), we deduce from (AC) and (A7) that 

a(u)b(v) = sinh(v — u) _1 (sinh(w — u + X)b(v)a(u) 

- e v ~ u sinh(A)6(«)o(«)) 2 , 

d(v)b(u) = y^ 2 smh(v — m) 1 ( sinh(w — it + X)b(u)d(v) 

- e u ~ v smh(X)b(v)d(u)) , 

where now 6(u) and commute. Defining the 'pseudo-vacuum' 



V« as 



.7=1 V ' 



the eigenvectors of T^^u) are given by 

F^' ni = b{u {2) ) . ..b{u%)F {im ® F<®. 



(A20) 
(A21) 

(A22) 
(A23) 
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FIG. 5. Action of T m (u) on F (1)(0) ® F (0) . 

Graphically, the action of T"W(u) on F^'' ni is depicted in Figure |j. This diagram arises from Figure^ when D(u) 
is dragged down through all B's and all unwanted diagrams are discarded. The edges on the diagonal on the left in 
Figure ||now represent the 'pseudo-vacuum' in the space Vnl ® ■ ■ • <8 Vj . The action of the bi on F^^ is 

represented by buil ding a lattic e on this diagonal in the same way as the operators Bi did on , see Figure § The 
commutation rules ( A20 ) and ( A21 ) can then be repr esent ed graphically similarly. 

The eigenvalues AW corresponding to the vectors (A23) are given by 



AW(tt) 



v 3 



'4yl) N yi 2n2 sinh(u) Ar J| sinh(u - u£> + A) x 



fc=i 



n 



(2) 

sinh(u z — u + A) 



• i./ ( 2 ) \ 
i=i svah(u) — u) 



+ v2 2ni {4yl) N yl ni - 2n2 



s\nh{u) N Y[ sinh(u — ujj. ) x 



fc=i 



n' 2 sinh(u — m[ 2 ' + A) 
i=i svah{u-u\') 



(A24) 



The eigenvalue (J,(U) is simply given by 



(A25) 



The 'unwanted terms' generated by the action of T^\u) on ( A23 ) can be re ad of f from ( |A20| ) and ( A21 ). Using 
the commutativity of b(u) and b(v) they can be shown to cancel and make (A23) an eigenvector precisely when 
A^(u^) = 0. The numbers {mP"*} therefore satisfy the equations 



n 



• \.f ( 2 ) Wl 



-j^ sinh(u 



(2) (2) 

3 U l 



A) 



k=i sinh(u} i; - u^> + A) " sinh(4 2) - uf> 



A) 



2 2\ N 

2 2 
V2 X 3 



2 2 



(A26) 



'sinh^W + A)' 



Knowing A^^u) the first set of equations as given by ( A16 ) becomes 

" fj£_Y 2 T[ s i uh{11 

sinh(u 



sinh(u 



2 2 \ N 



(i) 



A) 



\y 2 y 2 J k=1 sinh(it^ — u*jp 



"2 :„i,/„(l) 



A) 



n 



uf - A) 



sinh(u 



(i) 



(A27) 



The eigenvalue combining the expressions ( Alii ) and (A24) becomes 
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A(tt) = 2 " 1+2 " 2 ^ 2 " 2 sinh(« + A) W n 



sinh(u — u^ 1 ' — A) 



, ( 2 2\ N -2ni -2ri2 • \.r \N TT 

+ (^22/3) ys Vi smh(u) J { 



=1 sinh(u — u k 
sinh(u — + A) 



fc=l sinh(u-uj^) 



sinh(u — — A) 



^ sinlK// - ( /; J: i 



" 2 ' V, ( 2 ) _i_ \ \ 



1=1 smh(u-u\>) 



(A28) 
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